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densation in its vacuum. It consists of the color D4-branes and DO-branes as a background 
CO ■ 



We consider the model of holographic QCD with asymptotic freedom and gluon con- 



and the flavor D8-branes as a probe. By taking a specific field theory limit, the effective 
coupling decreases to vanish in UV region. We then introduce the uniformly distributed 
baryons in terms of the baryon vertices and study the density dependence of chiral con- 
densate, which is evaluated using the worldsheet instanton method. In the confined phase, 
the chiral condensate as a function of density monotonically decreases in high baryon den- 
sity. Such behavior is in agreement with the expectation, while in extremely low density it 
increases. We attribute this anomaly to the incorrect approximation of uniformity in very 
low density. In the deconfined phase the chiral condensate monotonically decreases in the 
whole region of density. 
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1. Introduction 

Chiral symmetry is one of the most important properties which control the behav- 
ior of the strongly interacting nucleon system and therefore studying its order parameter, 
the chiral condensate, especially in the dense matter is the most urgent and important 
to understand the properties of hadrons in nuclei, neutron stars and so on. So far many 
speculations have been made about the behavior of the chiral condensate in dense medium 
without definite proof. One of the intuitive and acceptable ones is that the chiral conden- 
sate has a finite value at zero density and vanishes at a certain high density where the 
chiral symmetry is restored, and that it interpolates between the two by monotonically 
decreasing behavior [1]. However, no one succeeded in proving this from a first principle 
in a theory with strong coupling since there has not been any reliable calculational tools 
for the strongly interacting system especially in the presence of the dense fermions. 

Since the gauge/gravity correspondence [2] is good for strongly interacting gauge 
theory, it is tantalizing to see what happen if we apply this method to QCD and such 
application is called holographic QCD. The flavor physics of QCD has been studied in 
many models of holographic QCD, for example, the D3/D7 model [3], the D4/D6 model 
[4] and the D4/D8 model [5,6]. One can incorporate baryons into those models by the use 
of the baryon vertices [7], the D5-branes wrapping S 5 in the D3-branes background or the 
D4-branes wrapping S 4 in the D4-branes background. Such models with baryonic medium 
have been attracting many interests [8-15]. However we have not found a model which has 
the behavior of chiral condensate in complete agreement with the speculation mentioned 
above. 

Although the D4/D8 model is one of the most successful models of holographic QCD, 
it has several shortcomings: current quark mass can not be incorporated, and as a con- 
sequence the chiral condensate can not be naturally introduced. Various prescriptions for 
this problem have been suggested by [16-22]. The method by [19] allows us to calculate 
the quark mass and chiral condensate (ipip) from the open Wilson line, whose expectation 
value is given by the minimal surface of the open string worldsheet surrounded by the 
D8-branes, the worldsheet instanton. Therefore the chiral condensate can be evaluated 
by the area of that minimal surface in the analogy of the computation of Wilson loop by 
[23]. Following this method, we have studied the chiral condensate in both confined and 
deconfined phases [15]. There, we found that the chiral condensate decreases in low density 
but increases in high density so that the low density behavior in this model agrees with 
the expectation, but the high density one does not, as in all other known models. 
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The D4/D6 and D4/D8 models use the common background yielded by the D4- 
branes, in which the gluon sector of four-dimensional dual gauge theory becomes non- 
supersymmetric due to an anti-periodic S l compactification of fermions [24]. This is an 
advantage of this background, because the ordinary QCD does not have supersymmetry. 
However the effective coupling e^ diverges as one goes to UV, namely, the theory has oppo- 
site behavior to asymptotic freedom. In this paper we improve this problem by introducing 
DO-branes smeared on the D4-branes. We take a specific field theory limit to achieve the 
asymptotic freedom. In general, the Dp-branes smeared on the D(p + 4)-branes play the 
role of instantons, so that the Dp-branes provide a gluon condensate to the field theory on 
D(p + 4)-branes [25]. 

The rest of the paper goes as follows. In Section 2 we suggest the new model whose 
background is given by taking the specific field theory limit of the supergravity solution of 
the D4-branes on which the DO-branes are smeared. In Section 3 we consider the D8-branes 
in the background of confined phase and then numerically calculate the chiral condensate 
depending on the density of baryons introduced as baryonic D4-branes (baryon vertices). 
In Section 4 we comment on the chiral condensate in deconfined phase which is analyzed 
in the similar way to the confined phase. Finally Section 5 is devoted to the conclusions. 

2. A new model of holographic QCD by D4-branes and DO-branes 

In this section, we propose a new model of holographic QCD which consists of N± 
color D4-branes, Nq DO-branes and Nf flavor D8 and anti-D8-branes. The DO-branes are 
smeared on the color D4-branes and induce a gluon condensate. The theory has U^N^) 
color gauge symmetry and U(Nf) x U(Nf) chiral symmetry. The flavor D-branes are 
treated as a probe embedded into the background produced by the color D-branes. 

2.1. D4+D0 solution in II A supergravity 

We consider JV4 D4-branes on which N$ (Euclidean) DO-branes are homogeneously 
smeared. The configuration of these D-branes are shown in Table 1. 



N4 D4-branes 
Nq DO-branes 



X° X 1 X 2 X s X 4 p 6 7 8 9 



Table 1 The configuration of the D4-branes and DO-branes. 
"~" stands for a smeared direction. 



In IIA supergravity the non-extremal black brane solution describing such D-branes has 
already been known [26-30]. The solution has the following metric and dilaton: 



ds 2 = ^dX,dX» + -^==(dX*) 2 + ^ThHl[^+p 2 d^ , (2.1a) 

e * = 9s {jh) 4 ' 9s:=e-^\ (2.16) 

where p = 0, 1, 2, 3 which are contracted by rj^ u = diag(— 1, 1, 1, 1). H$, H^ and F are the 
harmonic functions: 

#o(p) = l+(^) 3 , # 4 (p) = l+(^) 3 , F(p) = l-[P^\ (2.2) 

po and P4 are related with the horizon, p^, and the D-brane charges, Nq and N4, by 

Pi = -\pl + ^ + ^ 2 «(^^) 2 , Pi = ~\p\ + \j\pl + WsQlNl , (2.3) 

where V4 denotes the volume of the space on which the DO-branes are smeared, i.e., 
V 4 := JdX°dX 1 dX 2 dX 3 . 

There are two Ramond-Ramond (RR) fluxes, G^ and G(8)> which are associated with 
the D4-branes and DO-branes respectively. Since the fluxes obey (2tt£ s )~ 3 f Gu\ = N4 and 
(27i£ s )~ r f G(g) = iVo by the definition of the D-brane charges, one can write down the 
fluxes as 

= (27r4) 3 iV 4 rf04 ^ = (^sYNo dx0 a dxl a dx2 a dx3 a ^ ^ (2 4) 

O4 1/404 

where O4 is the volume of unit four-sphere, i.e., O4 = 87r 2 /3. 

2.2. The field theory limit 

Let us consider the Maldacena limit [31], 

■|r =: U is fixed, £ s ->■ . (2.5) 

Since t/ has the dimension of energy, the system under the limit (2.5) lies in the finite 
energy configuration. We additionally fix the ratio, 

(2n£ s ) 4 N 1 , N 
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to be finite. Since the leading terms of Hq and H 4 as £ s — > are 1 
the metric (2.1a) and the dilaton (2.16) are reduced to 



3 /7r^ s iV4 1 ( U h 
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(2.7) 



where £4, := Ph^ s 2 ■ Rescaling the radial coordinate back, we can rewrite (2.7) as 
ds 2 = V^dX,dX» + (|) 3 y^(^ 4 ) 2 + Mf){wp) + P 2dn " 



(2.8a) 
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e^ = a s ^(-J , F(p) = l-(y) , &:=itg a P a Nt. (2.86) 

Simultaneously the RR fluxes (2.4) are described as 



^^i d n t , G(s) = (?^i 3 ^-o.-. 



G(4) = if!!Mili dn4 , G{s) = \^J^^ dX ^dX l AdX 2 AdX 3 AdQ 4 . (2.9) 



We note that at the v — > limit the metric (2.8) does not become one in the field theory 
limit of N4 D4-branes (without a DO-brane). In order to obtain it, we have to set Nq = 
before taking the field theory limit. Therefore our D4+D0-branes background (2.8) is not 
valid at v = and is essentially different from the D4-branes background. 
For later convenience, we rescale the coordinates by R as 

a X ^ X 4 P 

R ' y R ' i?' 



If v is not fixed, the small £ s expansions of Hq and Ha become 

„ _ 1 1 K 2 g 2 J 2 s / (27r) 4 JV \ 2 4 -WsN±, n(f o, 

To the leading order of £ s , #0 is approximately equal to one, namely, the DO-branes do not 
affect the geometry, and the resulting metric and dilaton are nothing but those generated by only 
D4-branes. The metric containing the subleading term in Ho was studied by [32-34]. 



so that x M , y and r are dimensionless coordinates. Then the metric and dilaton (2.8) 
become 

R~ 2 ds 2 = ^dx.dx^ + r3fA L {r) dy 2 + ^ (^- + r 2 dnj) , (2.10a) 

V* 7 r \JM{r) J 

e * = 9s^/2> /M(r) = l-(-^J , ^ = irg a ^N A , (2.106) 



(27r£ s ) 4 iVo 1 . /"^ ^ 
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" = ^ 4 iV v*'-=Jdx= w , (2.10c) 

where tm := Ph/R, while the RR fluxes (2.9) are 

= (27r4) 3 iV 4 ^ 4 ^ ^ = (2ir£ s ) 3 R 4 N 4 u dx0 ^ ^ ^ ^ ^ ^ ^ ^ ^^ 



We compactify the y direction by a circle with a period /3 y in the same manner as [24] , 
so that we obtain the four-dimensional dual gauge theory. In order for the smoothness at 
t = tm, the period f3 y has to be 



fty 



°'M 



Then the Kaluza-Klein mass and the Yang-Mills coupling in the four-dimensional dual 
gauge theory are described as 

2n 3r 2 

Mkk = w = vfe ' g ™ = 2n9sisMKK • (2 - 12) 

The background given by (2.10) and (2.11) corresponds to the confined phase of the 
dual gauge theory. If we take the T-dual on the y direction, this D4+D0 becomes the 
D3+D(— 1) system which was studied by [25]. The D(— l)-brane plays a role of gluon 
condensation in the dual gauge theory, because a Dp-brane can be realized as an instanton 
on D(p + 4)-branes. Therefore the DO-branes smeared on the D4-branes behave as the 
instantons and provide the gluon condensation. 

We need to clarify the parameter region in which the supergravity approximation is 
valid. First, in order to suppress the loop correction, the effective string coupling must be 
small, namely, e^ <C 1. Since e^ decreases as r, this constraint is realized as o s z/ 3 ' 4 r M <C 
1, that is, 

v^A<iv|(4M K K) 2 , (2.13) 
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where A is the 't Hooft coupling, A := <7 Y m-^4- Second, the curvature must be smaller than 
£~ 2 for neglecting a higher curvature correction. We calculate the scalar curvature from 
the metric (2.10a), 

^) = ^2^(^M + r 3 ). (2.14) 

Since this is the increasing function for large r, it is necessary to introduce the UV cutoff 2 
Too so that 7^.(7-00) <C £~ 2 , namely, r^ <C -JvR 2 £~ 2 . This constraint with tm < r oo leads 
to 

(4M KK ) 2 «^A. (2.15) 

When we assume that 4-^kk is finite, (2.13) and (2.15) are consistent with the 't Hooft 
limit, i.e., the N 4 — > 00 limit with the large fixed A. Since we have taken the small £ s limit, 
the finiteness of £ s Mkk implies that Mkk is large, in other words, the compactification 
period (3 y is small. Therefore the dual gauge theory flows to be four-dimensional. We note 
that, in holographically computing physical values of the gauge theory, £ s and Mkk always 
appear in the combination of £ s Mkk- 

We can also write down the background geometry corresponding to the deconfined 
phase. By exchanging the coordinates and double Wick-rotation, namely, (x°, y) — > 
(iy,ix°), in (2.10) and (2.11), we obtain the metric, 

ir 2 ds 2 = _i?Mr), dx o)2 + ^t d$ 2 + dy 2x + V" (dr 2 + r2dQ2 \ 

\ru ' r 6 \f T {r) J 

0-13/2 > fT( r ) = 1 -(—) » R 3 = 7rg s £ 3 s N 4 , (2.166) 



R 4 v 4 N 4 
and the RR fluxes 



v 4 := I dx 1 dx' 2 dx 3 dy , (2.16c) 



(2^) 4 iVo 1 .. . / UjJjJ 



'1 



{2ti£ s ) 3 N 4 ir ^ „ {2tt£ s ) 3 R 4 N 4 u , , r2 ,3 , ^ , n ^ 

Gu) = - — -^ -dn 4 , G( 8 ) = - 1 - —dx 1 A dx 2 A dx 3 A dy A dQ 4 . (2.17) 

il 4 il 4 



In our field theory limit, the harmonic function Ho in (2.2) is reduced to Ho = vjr 3 . Let 
us naively imagine another limit where Hq has the form, Hq = 1 + vjr 3 (cf. [32,33]). The scalar 
curvature in this case is 

S3 ~6 



n{f) 



V0R 2 f 2 (l + r 3 /i>) 5 / 2 



:', -:■: r /~3 n ~3 s _^ ,, •:; c '.'. 

v 4j> 2 



M 



■if -2r M ) - -—(r -r M ) 



At small f, this curvature coincides with (2.14). In other words, our theory focuses on the small 
r region in the theory of Ho = 1 + vjr 3 . 
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The transition between the background geometry for the confined phase (2.10) and the 
one for the deconfined phase (2.16) is known as the Hawking-Page transition. Since we 
exchanged the time direction with the radial one, the configuration of D-branes is slightly 
modified as Table 2. 



N4 D4-branes 
Nq DO-branes 



„0 „,1 Jl _3 



x x x x y r 6 7 8 9 



o o o o o 

O r*-> r*-> <~*-' <~*-> 



Table 2 The configuration of the D4-branes and DO-branes for the deconfined phase. 

Since the geometry (2.16) is a black hole, we can introduce a temperature by compacti- 
fying the Euclidean time direction r (:= ix°) with a period f3 T which is determined to be 
(An/3)y/i'r^, 2 . Then it leads to the (Hawking) temperature, 

1 3r 2 

T= Wr = ^R- (2 ' 18) 

It is remarkable that the power of effective couplings e^, (2.106) and (2.166), has a 
negative value, —3/2. This implies the UV asymptotic freedom. On the other hand, the 
dilaton in the D4-branes background without DO-branes behaves like e^ ~ r + 3 / 4 [4-6] , that 
is, asymptotically non-free. Therefore we expect that the UV behavior in our holographic 
QCD model is improved by virtue of the DO-branes. 



3. Dense baryons in the confined phase 

3.1. The action of flavor D8-branes 

In order to consider the flavor dynamics in the confined phase, we embed Nf D8- 
branes and Nf anti-D8-branes into the background (2.10), in which the D8-branes and 
anti-D8-branes are connected with each other, so that they become Nf U-shape D8-branes 
(see fig. 1). This implies the chiral symmetry breaking from U(Nf) x U(Nf) to U(Nf). 



>/■ 




fig. 1 The embedding of D8-branes. 



The D8-branes are stretched on (x°, x 1 , x 2 , x 3 , r, S 4 ) and the collective coordinate y 
depends only on r (see Table 3). 
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Table 3 The embedding of the flavor D8-branes. 
Then the action of D8-branes without a worldvolume gauge field is 



S 8 = -T / d 4 xdrvJr 2 f M (r)(y'(r)) + 



v 



r A f M (r) 



(3.1) 



where T = Nffi8R 9 ^4/g s , A*8 = (2tt) 8 £ s 9 and ' denotes the derivative with respect to r. 
We calculate the equation of motion from the action (3.1), 



d 

d;r 



r 2 fM(r)y'(r)lr 2 f M (r)(y'(r)) + 



r A f M (r) 



-1/2' 



0. 



(3.2) 



Integrating this equation with respect to r and imposing the boundary condition, 



y(r ) = o, 2/(7-0) = 00 



(3.3) 



we obtain 



{y'(r)Y 



v 



r 2 f M {r) 



1 



(3.4) 



r 6 (/ M (r)) VTcJmM 

Therefore we can describe the separation between the D8-branes at the UV cutoff, r = r^, 
as 

(3.5) 



/ = 2/ dry'(r) = 2l ° dr -^—( ' C^X -\ 



to -/r r 3 f M (r) \r$f M (r ) 

Now let us turn on the baryons which are homogeneously distributed in M 3 (x 1 , x 2 , x 3 ). 
These baryons are realized by baryon vertices, that is, the (baryonic) D4-branes wrap- 
ping S 4 [7]. We assume that the baryon vertices are located on the D8-branes at 
r = r c . From the viewpoint of the worldvolume theory on the D8-branes, the number 
of baryons corresponds to U(l) charge, because the Chern-Simons action includes the 
term, A% (1) tr(F su ( N ^) 2 . 



More concretely we study the action of D8-branes. It consists of the DBI action and 
the Chern-Simons action, 

Sg = 5t)bi + Scs ■ 
The U(Nf) worldvolume gauge field, A, on the D8-branes can be decomposed as 

A = A su ^ + -^A u ^ . 

Since we are interested in the U(l) baryon charge, we focus on the time component of U(l) 
gauge field, A , which depends only on r, and we assume that the other components 
are equal to zero. For later convenience we define the dimensionless field by 

a(r) := ^L<« . (3.6) 

Due to this gauge field, the DBI action (3.1) is modified as 



// 9 1J I 9 

dxdr£ DB i, £DBi = v\r 2 f M (r)(y'(r)) + (a'(r)) • (3.7) 

y r 4 f M {r) r x 

We introduce the baryon vertices at r = r c . We then have to take into account the 
Chern-Simons term coupling with the four-form flux G(4), 

^ 8 (27T^ 



"I) 3 



3! /«(4)Ao; (5) = ^y'a; (5) , W(5) = Tr(AF 2 - \a*T + -^ 6 ) , (3-8) 
because this includes the term inducing the baryon charge, 

Scs = §^/ A<fr-|=^MF-^>) 2 . (3.9) 

Since the baryon vertices are at r = r c and are homogeneously distributed in R 3 , trF 2 is 
identified with the baryon number density as 

±_ tT ( F su(N f )f = nB6{r _ Tc) (3 1Q) 

Therefore (3.9) becomes 

S CS = T fd 4 xdr3N B a(r)5(r-r c ), N B :=n B {^-\ . (3.11) 

We note that, although there is the other RR flux, (?(§), in (2.11), it does not couple with 
a(r), the time component of U(l) gauge field. 3 



Compared with the source of baryon charge (3.8), (2n£ 2 ) 3 G(4) A A A T 2 ~ 0(£®), we can 
neglect (2tt^) 4 * G {8) A A^ {1) A T 3 ~ O^J 1 ) due to small 4- (See also (2.11).) 

9 



3.2. Force balance condition between D8-branes and baryonic D^-branes 

Since the baryonic D4-branes are attached on the D8-branes at r = r c , the force from 
D4-branes has to balance with that from the D8-branes [8]. 

We remind the total action of D8-branes given by (3.7) and (3.11 

S 8 = T d 4 xdr [-Xdbi + 3N B a(r)5(r - r c )] , 



i.e., 



C 



DBI 



v\ r 2 f M (r)(y'(r)) + 



v 



1 



r 4 fM(r) r 
The equations of motion for y(r) and a(r) are respectively 



--(a'(r))' 



(3.12a) 
(3.126) 



d_ 
dr 



r 2 fM(r)y'(r) 



0, 



£dbi 
d'{r) = 3N B S(r - r c ) , 

where we introduced the electric displacement field, 

5£ D bi v 2 d{r) 



d(r) :-- 



(3.13a) 
(3.136) 

(3.14) 



5a' {r) r£ D Bi 
From (3.136) we obtain 

d{r) = 3N B , 

that is, the electric displacement field becomes the constant proportional to the baryon 
density. With the boundary condition, y'(r Q ) = oo, we can solve (3.13a) as 

,2 v ( r 2 f M (r)(is 2 +rd 2 ) 



{y'(r)Y 



1 



(3.15) 



r 6 {f M (r)) 2 \rofM(r )(v 2 +r d 2 ) 
We are now considering the baryonic D4-branes which are attached on the D8-branes 
at r = r c . Then we can imagine two cases about the shape of D8-branes. These depend 
on the direction of force from the baryonic D4-branes. 






(a) (b) (c) 

fig. 2 The force balance between the D8-branes and baryonic D4-branes. 
(a) The V-shape D8-branes. (b,c) The W-shape D8-branes. 
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One is the case as fig. 2a where the V-shape D8-branes are pulled down by the falling 
D4-branes. We remember that this is the case for the generalized Sakai-Sugimoto model, 
i.e., the D4-branes background without the DO-branes [8]. The other is as fig. 2b where 
the W-shape D8-branes are hung up by the floating D4-branes. 

In order to clarify the shape of D8-branes in our D4+D0 model, we calculate the on- 
shell action of baryonic D4-branes. The baryon vertices with density ub are identified with 
the ubVs D4-branes wrapping S A . Note that v^ (:= f d 3 x) is the volume of R 3 . Therefore 
the on-shell action of baryonic D4-branes is 

5*4 = —(nBVs)/J,4 / dx°dfl4 e~^ a/— det #D4 =: — / dx° £4 , £4 = -Tv^d^ — . 

J on-shell J O \ r c 

(3.16) 
Since the energy, £4, decreases as r c , the D4-branes yield an upward force (to the positive 
direction of r). Therefore our D4+D0 model prefers the W-shape D8-brane than the 
V-shape one. We also need to take into account the possibility as fig. 2c, in which the 
W-shape D8-branes and the baryonic D4-branes are connected by fundamental strings 
and the position of D4-branes is rj, (> r c ) [35]. The action of the baryonic D4-branes and 
fundamental strings are described as 



S4+st = -(ubVs] 



/J4 / dx dQ 4 e ^ y/ - det gr>4 + --J2 \ dx ° / dr^- det String 



shell 



= -Tv 3 d^ I dx° —— + 



dr 



' /yt<J /Y' 



dx°£ 



4+st 



Since the derivative of the energy ^4+ s t with respect to r^ is always positive, ^4+ s t has 
a minimum at r\, = r c . In other words, the configuration like fig. 2b is favorable than 
that like fig. 2c. Finally the W-shape D8-branes on which the baryonic D4-branes directly 
attached (fig. 2c) is the (at least classically) stable configuration . 

For the W-shape D8-branes we calculate the separation at the UV cutoff r = r^ in 
terms of (3.15), 

l = 2(- p\ J rX \dr\y\r)\ 

r c cfoo\ rr. / Jit (Jl,.l 1 „j1\ \ -3 



r Jr 
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fig. 3 The force balance condition: /§ cos6> = f&. 

Now let us calculate the force balance condition between the D8-branes and the bary- 
onic D4-branes (see fig. 3). 4 The Legendre transform of the DBI action of D8-branes (3.7) 



is 



Sg := T dTxdr(a'(r)d + £ 



'DBI 



= T J d 4 xdrJ(v 2 + rd 2 )fr 2 f M (r)(y'(r)) 2 + r4/ ^ (r) ) ■ ( 3 - 18 ) 



From this equation, we can read that the effective tension at r = r c is 



Tv 3 i \v 2 

The angle between the D8-branes and the r axis at r = r c is computed as 

Jg^dr 



(3.19) 



cos6> 



^/g rr dr 2 + g yy dy* 



' 1 _ rlf M {r )(v 2 + r d 2 ) 
r 2 fM(r c )(v 2 +r c d 2 ) 



(3.20) 



The D8-branes provide the force, fs cos 9, along the r direction. On the other hand, the 
force from the D4-branes is evaluated in terms of (3.16), 



n 



1 d£ 4 



'Qrr 0!f c 



6R 



-Z/4 



z dy/f. 



M\Tc 



(3.21) 



As a result, from (3.19), (3.20) and (3.21), we obtain the force balance condition, /§ cos 9 
/ 4 , namely, 



y 



; 2 + r c d 2 
r c fM(r c ) 



L r 2 fM(ro)(v 2 + r d 2 ) \ 
V r 2 f M (r c )(u 2 +r c d 2 ) J ' 



(3.22) 



In this paper we calculate the forces in the left half (y(r) > 0). Therefore the actual total 
forces are twice and balance as 2f% cos 6 = If a- 
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3.3. Chiral condensate 

In order to consider the chiral condensate, we adopt the method suggested by [19], in 
which the open Wilson line operator, 

/1 11 

was introduced, i, j are indices of the fundamental representation of U(Nf). The AdS/CFT 
correspondence allows us to evaluate the vacuum expectation value of the open Wilson line 
operator by the use of the on-shell action of fundamental string whose boundary is the 
flavor D8-branes (see fig. 3), that is to say, 

(0>)~6{{0), (<D) = e- s °, S = ^Jj 2 o-V^~g, 

to leading order of a' (= £ 2 ). One can realize (O) as the chiral condensate {'ipip). So is, in 
other words, the area of the worldsheet instanton. By the use of (3.17), So can be written 
down as 



d2 fl/2 /'Too t->2 / rr c f roo \ 

So =<L *L* = i0.U n + Lr T -- TM (X23) 



We note that So does not have UV divergence, because we have introduced the UV cutoff, 
Too, in our model as we have discussed in Section 2. However So still has a large value 
depending on the cutoff. So we consider (O) divided by the large constant, 



(C)oo := exp 



R 



2 



27r£? 



l(r oo - r M ) 



(3.24) 



Let us evaluate the chiral condensate numerically. In order to sweep out tm from 
numerical computations, we rescale the variables as 



w :— — , d:=JrMol, I :— r 2 M l 



rM 



Then the force balance condition (3.22) is rewritten as 

Wq(1 - Wq 3 )(u 2 + w d 2 ) = w 2 (l - w~ 3 ) ( v 2 + w c d 2 -^— ) =: Q(v,d,w c 
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This equation determines wq as a function of u, d and w c . From (3.17), (3.23) and (3.24), 
the separation and chiral condensate are described as 



, = 21 f + fk-#-f^^±^-l 
I wo J wo J w s -l\ Q(u,d,w c ) 

( R 2 



(3.25) 



(g) 

(^)oc 

J = 



exp 



W^M , 



W c fW L 

+ / w 

WQ J WQ 



v fw 2 (l-w s )(v 2 + wd 2 ) 



W 2 + W + 1 



Q(u,d,w c 



1 . (3.26) 



By fixing u, I and w^, we calculate the density d dependence of exp(J). 
exp(7) exp(J) 





exp(I) 



exp(J) 




v = 10 



v = 1 



i/ = 0.1 



rd 



(c) (d) 

fig. 4 The plots of the density d dependence of exp(J) with w^ = 10 4 and / = 0.1 fixed, 
(a) v = 0.1. (b) v = 1. (c) i/ = 10. 

(d) The plots for ^=0.1 (solid), 1 (dotted), 10 (dashed). Every exp(J) converges to 
a certain constant at d — > oo. 



For instance, when we set w^ = 10 4 , / = 0.1 and v = 0.1, 1, 10, we obtain fig. 4. Note 
that the contribution from large w to the integrations in (3.25) and (3.26) is very small, 
therefore exp(J) is not sensitive to the cutoff w^. We can see in fig. 4 that the chiral 
condensate, (O) (= (0) QO (e I ) R /( ni s r M)^ increases in very low density. This is different 
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from our intuition from ordinary QCD. On the other hand, in high density, the chiral 
condensate is a monotonically decreasing function of d and converges to a finite value as 
d — > oo. In all other known models of holographic QCD [12-15], the chiral condensate 
increases in high density and it is completely opposite to the expectation from ordinary 
QCD. Therefore one can say that our model drastically improves the high density behavior 
of chiral condensate. 

4. Comments on the deconfined phase 

4-1. The action of flavor D8-branes 

The background corresponding to the deconfined phase has been given by (2.16) and 
(2.17). We embed the Nf D8-branes (and Nf anti-D8-branes) into it in the same manner 
as we did in the confined phase, that is to say, we identify the worldvolume coordinates 
of the D8-branes with (a; ,^ 1 , x 2 , x 3 ,r, 4 ) and assume that the collective coordinate, y, 
depends only on r. Since we are interested in the U(l) baryon charge, we take into account 
the U(l) gauge field denoted by (3.6). Then the DBI action is described as 



4 d Bi° = -Tj#xdrrtgg , 4 d Bi } = ^rif T (r)(y>(r)) 2 - 1 -(a>(r)) 2 + 1 -. (4.1) 

There are two Chern-Simons terms coupling with the U(l) gauge field, a(r). One is the 
coupling with the four-form flux (2.17), i.e., f G^\ Aw( 5 ), originated from the background 
D4-branes. As we explained in Section 3.1, U($) includes A u ^ 1 'tr(F su ^ Nf ') 2 . Therefore it 
provides us the source term of the U(l) baryon charge, 

S@ =T f d 4 xdr 3N B a(r)5(r - r c ) , (4.2) 

in the same way as (3.11). This also implies that the Ng baryonic D4-branes wrapping S 4 
are located at r = r c . 

In the deconfined phase there is the other Chern-Simons term induced by the eight- 
form flux (2.17), G(8), of the DO-branes, namely, f G^g) /\Aodx°. This term includes 



q (D0) 



r/A*W(rWO. (4-3) 



We note that this Chern-Simons term yields the U(l) charge even if there is not a baryon 
(N B = 0). 
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4-2. Without a baryon 

In this subsection we concentrate on the case of no baryon, that is, Nb = 0. From 
(4.1) and (4.3), the action of D8-branes is 



S. 



(dec) I 



I7V B =0 



-T / <i 4 £<ir [.Cqbi — 3vy'(r)a(r)] . 



We calculate the equations of motion from this action, 



d 

dr 

d 
dr 



ur 2 f T (r)y'(r) 



C 



(dec) 



— 3a(r) 



DBI 



va!{r) 



.rC 



(dec) 
DBI 



3y(r) 



0. 



(4.4) 

(4.5a) 
(4.56) 



Since these equations are symmetric under the constant shifts of the fields, we redefined 
the fields by 

<i/,(r) := y(r) + d , a(r) := a{r) + C 2 . (4.6) 

Then, integrating (4.5) with respect to r, we obtain 

3a(r) 



ifj'(r) 



a'(r) 



Jr 3 /T(r){9r3/ T (r)(^(r)) 2 - 9(a(r)) 2 + r 2 f T (r] 

3\/r 3 / T (r)^(r) 

^9r 3 / T (r)(^(r)) 2 -9(«(r)) 2 + r 2 / T (r)' 



(4.7a) 



(4.76) 



where the integration constants are absorbed by C\ and C 2 - Note that, when (ip(r), a(r)) 
is a solution of (4.7), (— i()(r), —a(r)) is also a solution. 

The trivial solution is a(r) = ip(r) = 0. When we put a pair of boundary conditions, 
y( r oo) = =W/2 5 this solution leads to a pair of solutions, y(r) = ±//2, which represents the 
parallel D8-branes and anti-D8-branes separated from each other by /. This configuration 
preserves the chiral symmetry, U(Nf) x U(Nf). 

Let us look for a non-trivial solution. In order to find a U-shape solution whose tip is 
located at r = ro, we impose the boundary condition: 



V>(r o ) = 0, i/)'(r ) = oo 



(4.8) 



From (4.7a), one can interpret ^'{tq) = oo to r^f^o) — 9(a(ro)) = 0. Therefore the 
boundary condition (4.8) can be replaced with 



ip(r o ) = 0, a(r ) = -roVf( r o)- 



(4.9) 
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We are now considering the Nb = case, i.e., there is no source induced by a baryon. 
Therefore there does not exit a cusp singularity on the D8-branes. In terms of the solution 
of (4.7) and (4.9), such a smooth U-shape D8-brane is described as 



y {r) — ±ip(r) , a(r) = ±(a(r) - a(r )) 



(4.10) 



where C\ and Ci in (4.6) are determined so that y(r) and a(r(y)) are smooth at y(r ) = 0. 



y 

l.n - 



y 



40 -6a 



100 120 ' 



(a) (b) 

fig. 5 The plots of y(r) and a{r) with fixed Tt — 1 for ro= 2 (solid), 10 (dotted) and 
50 (dashed), (a) The plots of y(r). (b) The parametric plots of (a^.y^)). 

Since it is difficult to solve (4.7) analytically, we perform numerical computation. For 
instance, in fixing ty = 1, we obtain the solutions depicted in fig. 5. The D8-branes and 
anti-D8-branes are smoothly connected at r = ro, so that the U-shape D8-branes appear. 
This implies the chiral symmetry breaking from U(Nf) x U(Nf) to U(Nf). Although 
the gauge field diverges at r — > oo, we neglect the large r region. Because, as we dis- 
cussed in Section 2.2, the UV cutoff r^ is necessary due to the validity of supergravity 
approximation. 



4-3. With baryons 

We shall turn on the baryons which are the ubV3 baryonic D4-branes located at 
r = r c . The total action of D8-branes is given by (4.1), (4.2) and (4.3) because of iVg ^ 0. 
Although the equation of motion for y(r) is same as (4.5a), the one for a{r) is slightly 
different from (4.56) due to the Chern-Simons term (4.2). By redefining 



d(r) :-- 



5C 



(dec) 
DBI 



8a' (r) 
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- 3vy(r) 



(4.11) 



the equation of motion for a(r) has a simple expression, d'(r) = 3NsS(r — r c ). Therefore 
d(r) is a constant equal to 3Nb- 

In terms of (2.16) we write down the on-shell action of the baryonic D4-branes wrap- 
ping 5 4 as 



S. 



(dec) 



M<ll,t ' ) ' ( ' M = iTvsdrcy/frfc 



~ i dx"£r" } , £4 



(4.12) 



.(dec) 



Since the energy E\ monotonically increases as r c , the baryonic D4-branes obtain a 
force toward the negative direction of r. Therefore this force from the D4-branes balances 
with the force from the V-shape D8-branes in a manner as fig. 2a. By using the solution 
ofip(r) and a(r) that we obtained in the previous subsection, we can describe the V-shape 
D8-brane as 

y(r) = ±(ip(r) - ip(r c )) , a(r) = ±(a(r) - a(r )) . (4.13) 

Note that the gauge field a(r) is not continuous at r = r c because there is the baryonic 
D4-branes, namely, the source of U(l) charge. 

The Legendre transformed action of D8-branes is 



^dec) =? - f d 4 xdr ( u 2 + r ( rf + 3uy(r)) 2 )(r 2 f T ( r )(y'(r)) 2 + -} , (4.14) 

up to total derivative. From (4.12) and (4.14) we calculate the force balance condition, 



-dr c 
3 



2 Vr c 



/r(r c )(z/ 2 + r c d 2 ) 



rlfT{r c )(y'(r c )) + r c X 
The separation of V-shape D8-branes at r = r^ is 

1 = 2 I | 2 / / (r)| = 2(^(r 0O )-^(r c )). 



(4.15) 



(4.16) 



In the same way as Section 3.3, we shall study the density dependence of the chiral con- 
densate (O) by 



(g) 

(O)oc 



exp 



exp 



*£* 



R 2 ^ 

2n£ 2 s J 



^oo ry(r) 

dr J dy 
Jo 



v 



B? 



1(D) 



l (D) 



r 3 fT(r) 7t£ 2 

Mr) -I 



roc rl/2 

dr I dy 
Jo 



v 



r 3 f T (r) 



dr- 



(4.17) 



T u — Tn 
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exp(I(£>)) 




-iad 



fig. 6 The plot of exp(i7m) depending on the density d with fixed / = 0.5 
and ry = v = 1. 

For instance we numerically evaluate exp(/(£))) with fixed / = 0.5 and r? — v — 1, 
which is depicted in fig. 6. It shows that the chiral condensate monotonically decreases as 
the density d. This agrees very well with our expectation from the ordinary QCD. The 
exp(/(£))) for / = 0.5 and ty = v = 1 is approximately equal to 0.966 at d — 0, and 
converges to a finite value, 0.947, at d = oo (see fig. 6). 

5. Conclusions 

In this paper we proposed a model of holographic QCD whose background has the 
asymptotic freedom. It is the specific field theory limit of the supergravity solution of the 
color D4-branes on which the DO-branes are smeared. In taking the limit, we fixed the 
ratio v to be finite. Into this background we have introduced the flavor D8-branes and 
the baryonic D4-branes in order to study the chiral condensate in dense baryonic medium. 
The dilaton in the background of D4-branes without DO-branes behaves like e^ ~ r 3 / 4 . 
On the other hand, the dilaton in our background behaves like e^ ~ r - 3 / 2 5 that is, the 
effective coupling is asymptotically free. Although it is still different from (logr) -1 in 
ordinary QCD, the DO-branes improve the asymptotic freedom in our model. 

We have studied the chiral condensate by use of the open Wilson line, which we can 
evaluate by the area of open string worldsheet whose boundary is the D8-branes. In the 
confined phase the baryonic D4-branes balance with the D8-branes which have W-shape 
rather than V-shape. This is because the force from the D4-branes is toward the positive 
direction of the radial coordinate. This can be understood that the curvature increases for 
large radius so that the gravity is acting on baryon vertex upward. 

Under this force balance condition we have calculated the chiral condensate concretely 
with the fixed /, the separation between the D8-branes at the UV cutoff. The results for 
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three cases of v are shown in fig. 4. All of the three have the common features as follows. 
In the region of very low density of baryons, the chiral condensate increases as the density. 
This behavior is not what one expects from the QCD. We attribute this anomaly to the 
approximation of the uniform distribution of baryons in the space M 3 , which is not a good 
one when the density of baryon is small. In high density we have obtained a remarkable 
result, that is, the chiral condensate monotonically decreases. This is in good agreement 
with our intuition from the ordinary QCD. When the density goes to infinity, the chiral 
condensate converges to a certain constant. 

We have studied also the deconfined phase. We numerically solved the equations of 
motion without the baryons, and obtained the U-shape solution of D8-branes. It seems 
that the solution of the gauge field has UV divergence. However we do not need to seriously 
worry about the divergence, because basically the UV cutoff is necessary in our background 
due to the validity of supergravity approximation of bulk theory. Turning on baryons, we 
have calculated the chiral condensate. It is different from the confined phase that the 
D8-brane balancing with the baryonic D4-branes has V-shape in the deconfined phase. As 
a result, it monotonically decreases in whole range of the density (see fig. 6). This is in 
agreement with the expectation from ordinary QCD. In the deconfined phase there is the 
solution of parallel D8-branes, which preserves the chiral symmetry. Therefore we expect 
the chiral symmetry restoration occurs in certain high density [36]. One can guess that 
this transition is clarified by comparing the free energies of two configurations, namely, 
the V-shape D8-branes with baryon vertices and the parallel D8-branes. However there is 
subtlety due to the UV cutoff. Therefore we need to develop a way beyond the UV cutoff. 



Acknowledgement 

This work was supported by the NRF grant through the SRC program CQUeST 
with grant number 2005-0049409 and partially supported by the WCU project of Korean 
Ministry of Education, Science and Technology (R33-2008-000-10087-0). 



20 



References 

[1] G. E. Brown and M. Rho, "Scaling effective Lagrangians in a dense medium," Phys. 

Rev. Lett. 66 (1991) 2720. 
[2] J. M. Maldacena, "The large N limit of superconformal field theories and supergrav- 

ity," Adv. Theor. Math. Phys. 2 (1998) 231. [hep-th/9711200]. 
[3] A. Karch and E. Katz, "Adding flavor to AdS/CFT," JHEP 0206 (2002) 043. [hep- 

th/0205236]. 
[4] M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, "Towards a holographic 

dual of large N c QCD," JHEP 0405 (2004) 041. [hep-th/0311270]. 
[5] T. Sakai and S. Sugimoto, "Low energy hadron physics in holographic QCD," Prog. 

Theor. Phys. 113 (2005) 843. [hep-th/0412141]. 
[6] T. Sakai and S. Sugimoto, "More on a holographic dual of QCD," Prog. Theor. Phys. 

114 (2005) 1083. [hep-th/0507073]. 
[7] E. Witten, "Baryons and branes in anti-de Sitter space," JHEP 9807 (1998) 006. 

[hep-th/9805112]. 
[8] O. Bergman, G. Lifschytz and M. Lippert, "Holographic nuclear physics," JHEP 0711 

(2007) 056. [arXiv:0708.0326 [hep-th]]. 

[9] K. -Y. Kim, S. -J. Sin and I. Zahed, "Dense hadronic matter in holographic QCD," 

[hep-th/0608046]. 
[10] K. -Y. Kim, S. -J. Sin and I. Zahed, "The chiral model of Sakai-Sugimoto at finite 

baryon density," JHEP 0801 (2008) 002. [arXiv:0708.1469 [hep-th]]. 
[11] K. -Y. Kim, S. -J. Sin and I. Zahed, "Dense holographic QCD in the Wigner-Seitz 

approximation," JHEP 0809 (2008) 001. [arXiv:0712.1582 [hep-th]]. 
[12] S. Nakamura, Y. Seo, S. -J. Sin and K. P. Yogendran, "A new phase at finite quark 

density from AdS/CFT," J. Korean Phys. Soc. 52 (2008) 1734. [hep-th/0611021]. 
[13] Y. Seo and S. -J. Sin, "Baryon mass in medium with holographic QCD," JHEP 0804 

(2008) 010. [arXiv:0802.0568 [hep-th]]. 

[14] M. Edalati and J. F. Vazquez-Poritz, "Chiral condensates in finite density holographic 

NJL model from string worldsheets," [arXiv:0906.5336 [hep-th]]. 
[15] S. Seki and S. -J. Sin, "Chiral condensate in holographic QCD with baryon density," 

JHEP 1208 (2012) 009. [arXiv:1206.5897 [hep-th]]. 
[16] R. Casero, E. Kiritsis and A. Paredes, "Chiral symmetry breaking as open string 

tachyon condensation," Nucl. Phys. B 787 (2007) 98. [hep-th/0702155 [hep-th]]. 
[17] K. Hashimoto, T. Hirayama and A. Miwa, "Holographic QCD and pion mass," JHEP 

0706 (2007) 020. [hep-th/0703024 [hep-th]]. 
[18] O. Bergman, S. Seki and J. Sonnenschein, "Quark mass and condensate in HQCD," 

JHEP 0712 (2007) 037. [arXiv:0708.2839 [hep-th]]. 

21 



[19] O. Aharony and D. Kutasov, "Holographic duals of long open strings," Phys. Rev. D 

78 (2008) 026005. [arXiv:0803.3547 [hep-th]]. 
[20] K. Hashimoto, T. Hirayama, F. -L. Lin and H. -U. Yee, "Quark mass deformation of 

holographic massless QCD," JHEP 0807 (2008) 089. [arXiv:0803.4192 [hep-th]]. 
[21] A. Dhar and P. Nag, "Tachyon condensation and quark mass in modified Sakai- 

Sugimoto model," Phys. Rev. D 78 (2008) 066021. [arXiv:0804.4807 [hep-th]]. 
[22] S. Seki, "Intersecting D4-branes model of holographic QCD and tachyon condensa- 
tion," JHEP 1007 (2010) 091. [arXiv:1003.2971 [hep-th]]. 
[23] J. M. Maldacena, "Wilson loops in large N field theories," Phys. Rev. Lett. 80 (1998) 

4859. [hep-th/9803002]. 
[24] E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in gauge 

theories," Adv. Theor. Math. Phys. 2 (1998) 505. [hep-th/9803131]. 
[25] H. Liu and A. A. Tseytlin, "D3-brane-D-instanton configuration and M = 4 su- 

perYM theory in constant selfdual background," Nucl. Phys. B 553 (1999) 231. [hep- 

th/9903091]. 
[26] M. Cvetic and A. A. Tseytlin, "Nonextreme black holes from nonextreme intersecting 

M-branes," Nucl. Phys. B 478 (1996) 181. [hep-th/9606033]. 
[27] M. S. Costa, "Composite M-branes," Nucl. Phys. B 490 (1997) 202. [hep-th/9609181]. 
[28] I. Y. Aref'eva, K. S. Viswanathan and I. V. Volovich, "p-brane solutions in diverse 

dimensions," Phys. Rev. D 55 (1997) 4748. [hep-th/9609225]. 
[29] M. S. Costa, "Black composite M-branes," Nucl. Phys. B 495 (1997) 195. [hep- 

th/9610138]. 
[30] I. Aref'eva and A. Volovich, "Composite p-branes in diverse dimensions," Class. 

Quant. Grav. 14 (1997) 2991. [hep-th/9611026]. 
[31] N. Itzhaki, J. M. Maldacena, J. Sonnenschein and S. Yankielowicz, "Supergravity and 

the large N limit of theories with sixteen supercharges," Phys. Rev. D 58 (1998) 

046004. [hep-th/9802042]. 
[32] J. L. F. Barbon and A. Pasquinucci, "Aspects of instanton dynamics in AdS/CFT 

duality," Phys. Lett. B 458 (1999) 288. [hep-th/9904190]. 
[33] K. Suzuki, "D0-D4 system and QCD 3+ i," Phys. Rev. D 63 (2001) 084011. [hep- 

th/0001057]. 
[34] C. Wu, Z. Xiao and D. Zhou, "Sakai-Sugimoto model in D0-D4 background," 

[arXiv:1304.2111 [hep-th]]. 
[35] S. Seki and J. Sonnenschein, "Comments on baryons in holographic QCD," JHEP 

0901 (2009) 053. [arXiv:0810.1633 [hep-th]]. 
[36] O. Aharony, J. Sonnenschein and S. Yankielowicz, "A holographic model of decon- 

finement and chiral symmetry restoration," Annals Phys. 322 (2007) 1420. [hep- 

th/0604161]. 



22 



